Specific Heat and Tiiermodynamic Critical Field for Calcium under the Pressure at 

120 GPa 
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The free energy difference between the superconducting and normal state for Calcium under the 
pressure at 120 GPa has been determined. The numerical calculations have been made in the 
framework of the imaginary axis Eliashberg approach. On the basis of the obtained results the 
specific heat in the superconducting C'^ (T) and normal (T) state, as well as, the thermody- 
namic critical field He (T) have been obtained. It has been shown that the characteristic values 
of the considered thermodynamic quantities do not obey the BCS universal laws. In particular, 
AC (Tc) IC^ {Tc) = 2.48 and TcC^ (To) /H^ (0) = 0.154. 

PACS numbers: 74.20. Fg, 74.25.Bt, 74.25. Ha 



I. INTRODUCTION 

The superconducting properties of Calcium under the 
high pressure (») are intensively studied since 1981 (Dunn 
and Bundy) [l[. The pressure dependence of the critical 
temperature (Tc) has been determined by Yabuuchi et 
al. in 2006 [4]. The obtained results showed that Tc 
increases significantly with the increasing pressure from 
3 K to 23 K for p e (58, 113) GPa. Above 113 GPa the 
grow of Tc is considerably slower and at 161 GPa the 
critical temperature reaches the maximum value, which 
is equal to 25 K. We notice that in the considered pres- 
sure region. Calcium shows the complicated structural 
phase transitions [ll-[5|. The proposed structural phase 
diagrams for Calcium the reader can find in @ . 

In the presented paper we have calculated the free 
energy difference between the superconducting and nor- 
mal state for Calcium under the pressure at 120 GPa 
{Tc = 24 K). Next, the specific heat and the thermody- 
namic critical field have been determined. The numerical 
analysis was based on the Eliashberg equations on the 
imaginary axis [7]. 

Let us pay attention that the Eliashberg approach ex- 
tends the original idea of Bardeen, Cooper and Schrieffer 
, taking exactly into consideration the electron-phonon 
interaction. In the framework of the Eliashberg formal- 
ism, the strong coupling corrections to the BCS results 
are dependent on the value of the parameter ksTcluJin- 
The symbol win is called the logarithmic phonon fre- 



quency and Win — exp ^ 



In (17) 



For 



Calcium, the Eliashberg function (a^F(ri)) has been cal- 
culated in the paper Q , the maximum phonon frequency 
(i^max) and the electron-phonon coupling constant (A) 
are equal to 61.68 meV and 1.3 respectively. In the case 
of the BCS limit, the Eliashberg function is non-zero only 
for very high frequency, so that /csTc/win — > 0. In Cal- 



cium, value of the ratio kBTc/uJin is equal to 0.082. In 
this case the thermodynamic parameters can't be calcu- 
lated exactly in the framework of the BCS model. 



II. THE ELIASHBERG EQUATIONS 

The Eliashberg equations on the imaginary axis can be 
written in the following form [7] : 



M 



and 



K{n,Tn) - ii*e{ujc~ |w„J) 



A„, (1) 
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K (n, m) 



(2) 



where the symbol A„ = A (jaj„) denotes the order pa- 
rameter and Zn = Z {iuJn) is the wave function renor- 
malization factor; n-th Matsubara frequency is defined 
as: Wn = ^ (2n — 1), where (3 = 1/kBT. The electron- 
phonon pairing kernel K (n, m) is given by: 



K (n, m) 



dQ- 



-a'F{n). (3) 



In Fig. [T] we have presented the form of K (n, m) for the 
positive Matsubara frequencies and the temperature 1.16 
K. It is easy to notice that the pairing kernel is always 
positive and it achieves the strong maximum for w„ — 
LOm- The above result means, that Eqs. ^ and ^ can 
have the superconducting solution (A„ 7^ 0). 

In the framework of the Eliashberg formalism, 
depairing electronic interaction is described by 



the 
the 
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Coulomb pseudopotential: /i* = /i 

where ^ is defined by: ^ = ^(O)^^- The symbol p (0) 
denotes the value of the electronic density of states at 
the Fermi energy and Vc is the Coulomb potential. The 
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FIG. 1: The pairing kernel K (n, m) as a function of the num- 
bers n and m. 



quantities up and lod are the electronic plasma frequency 
and the Debye phonon frequency respectively - We have 
calculated the value of the Coulomb pseudopotential in 
the paper The following result has been obtained: 

/i* = 0.215. Finally, Q is the Heviside unit function and 
iOc represents the cut-off frequency: lOc = Silmax- 

The Eliashberg equations have been solved for 2201 
Matsubara frequencies {M = 1100) by using the method 
presented in [llj and [l^]. In the considered case, the 
obtained Eliashberg solutions are stable for T > 1.16 K. 



III. THE NUMERICAL RESULTS 

The free energy difference between the superconduct- 
ing and normal state (Ai^) for an interacting electron- 
phonon systems should be determined by using the ex- 
pression |13j : 



W) 



zvr X ^ 



n=l 



A2 



(4) 



where Zf^ and denote the wave function renormaliza- 
tion factors for the superconducting (5*) and normal (N) 
state respectively. In Fig. [5] we have plotted the depen- 
dence of (Ai^) on the temperature. From the physical 
point of view, the negative values of AF prove that the 
superconducting state is stable to the critical tempera- 
ture. 

We notice that the knowledge of the function Z^ and 
A„ enables also the determination of the open form of 
the electronic self energy for the superconducting state 
{T,^ (n)) 0. In the framework of the Eliashberg formal- 
ism, T,^ {n) is represented by the 2x2 matrix. For the 
half-filled energy band, the diagonal elements of Y,^ (n) 
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FIG. 2: The free energy difference between the superconduct- 
ing and normal state as a function of the temperature. The 
insets show the diagonal and non-diagonal elements of the 
matrix self energy on the imaginary axis. 



are imaginary and the non-diagonal elements are real. 
Additionally, it is occurring: Ef^ (n) = {n) and 
(") = ^21 ("•) • The dependence of Im [T,f^ (n)] and 
Re [Ef2 C*^)] on the positive values of the number n is 
presented in the insets in Fig. [21 It is easy to see that the 
values of Im[I]f]^ (n)] can be relatively high in compar- 
ison with Re[Ef2 (^)] ^^^d the considered function does 
not saturate for the large values of n. In the second case, 
we have found that Re[Sf2 (n)] has approximately a form 
of the Lorenz function. 
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FIG. 3: The specific heat for the superconducting and normal 
state as a function of the temperature. 

The difference in the specific heat between the super- 
conducting and normal state (AC = — C^) should 



be calculated by using the expression: ^C!{T) _ 

1 d^[AF/p{0)] 



P d{kBT) 



kBPlO) 

On the other hand, the specific heat in 
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FIG. 4: The thermodynamic critical field as ; 
temperature. 



function of the 



the normal state can be obtained with the help of the 



formula: 



C" (T) 
fcBP(O) 



is given by: 7 = |7r 



, where the Sommerfeld constant 
(1 + A). The dependencies of the 



specific heats on the temperature have been presented 
in Fig. [3l It is easy to see that at low temperatures, 
the specific heat in the superconducting state is expo- 
nentially suppressed. For the higher temperatures 
rapidly increases and the values of the specific heat in the 
superconducting state much exceed the values of C'^ . At 
the critical temperature the characteristic jump has been 
marked by the vertical line. 

Below we have calculated the values of the ther- 
modynamic critical field (cgs units): — 



temperature dependence 



of 



v/-87r [AF/p (0)]. The 

He I \/ P (0) has been shown in Fig. 21 

The knowledge of the thermodynamic functions C"^, 



and He enables the determination of the fundamen- 
tal ratios: 



and 



r2 



c) 



(5) 



where He (0) ~ He [T ^ l.\&K). For Calcium under 
the pressure at 120 GPa we have obtained: ri = 2.48 
and r2 = 0.154. Let us notice that the above results 
strongly differ from the canonical BCS predictions. In 
particular, [rijg^g = 1.43 and [r2]Bcg = 0.168 ^8,]. 



IV. CONCLUDING REMARKS 

The imaginary axis Eliashberg equations for Calcium 
under the pressure at 120 GPa have been exactly solved 
in the paper. On the basis of the obtained results we 
have calculated the free energy difference between the su- 
perconducting and normal state by using the expression 
given by Bardeen and Stephen. Next, the specific heat 
for the superconducting and normal state, as well as, the 
critical field have been determined. It has been shown 
that the ratios between the characteristic values of the 
calculated thermodynamic functions strongly differ from 
the values predicted by the BCS model. In particular, 
AC (Tc) /C^ (Tc) - 2.48 and TcC^ {Tc) I Hi (0) = 
0.154. 
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